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Abstract 

Within the framework of Weyl calculus we establish a quantum-classical corre- 
spondence for the time evolution of observables generated by a Dirac-Hamiltonian. 
This includes a semiclassical separation of particles and antiparticles. We then prove 
quantum ergodicity for Dirac-Hamiltonians under the condition that a skew prod- 
uct of the classical relativistic translational motion and relativistic spin precession is 
ergodic. 
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1 Introduction 



Dynamical properties of quantum systems have recently attracted considerable attention, 
in particular in connection with the question of quantum chaos. In this context a central 
problem is to characterise eigenstates of a quantum Hamiltonian semiclassically in terms of 
the dynamical behaviour of the associated classical system. Phase space lifts of eigenfunc- 
tions have proved to be especially suited for such investigations since in the semiclassical 
limit they converge to classically invariant measures on phase space. This property follows 
from a dynamical version of the correspondence principle, which is established in a math- 
ematically rigorous form through Egorov's theorem ||. Prominent examples of invariant 
measures on phase space are ergodic ones. On the quantum side a classically ergodic Li- 
ouville (i.e. microcanonical) measure corresponds to the property of quantum ergodicity. 
This concept goes back to Shnirelman |T^] and denotes the semiclassical convergence of 



the phase space lifts of almost all quantum eigenfunctions to Liouville measure; it has been 
proven in several situations [| [Tcfl . 

Quantum systems that possess spin degrees of freedom in addition to their translational 
ones require an extension of Egorov's theorem and of the concept of quantum ergodicity. 
In previous studies of non-relativistic quantum systems with spin [|T], |3|, whose dynamics 
are generated by a Pauli-type operator, we noticed that in general classical ergodicity 
of the translational motion is not sufficient to guarantee quantum ergodicity; the spin 
dynamics have to be considered as well. Here we present an extension of these investigations 
to Dirac-Hamiltonians. The main difference to the non-relativistic situation lies in the 
coexistence of particle and antiparticle states. This requires a separation of the Hilbert 
space into two subspaces that are almost invariant under the quantum time evolution. 
Only in these subspaces the usual semiclassical methods can be applied. As a consequence, 
quantum ergodicity for Dirac-Hamiltonians is concerned with projections of eigenspinors 
to the almost invariant subspaces, which in general are only approximate eigenspinors of 
the Hamiltonian. 

Due to lack of space we cannot present proofs here. Together with further references 
these can be found in , where we have investigated quantum ergodicity for more general 
matrix valued Hamiltonians. 

2 Background 

We consider the Dirac equation 

iH^(t,x)=H D if>(jt,x) (2.1) 
at 

in a fixed inertial frame in which the potentials A and <fi are time-independent, 

H d = col- (yV - ~ c M x )) + mc 2 p + e(f){x) , (2.2) 
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and the Clifford algebra, realised in the usual Dirac representation, is generated by a and 
(3. The propagation of Dirac spinors prescribed by (]2.1|) therefore takes place in the Hilbert 
space H = L 2 (R 3 ) <g> C 4 . 

For the purpose of semiclassical approximations it is useful to represent quantum ob- 



servables as Weyl operators. In the case of the Dirac-Hamiltonian ( |2.2| ) this means 

(iW)0*0 = {ov w [H D ty)(x) = j^— J j eiP-i-y) H D (p,^) dy&p . (2.3) 

Here the Weyl symbol H D (p,x) G C°°(R 3 x R 3 ) ® C 4x4 is the smooth function on the 
classical phase space R 3 x R 3 that arises from ( |2.2| ) by replacing j V with p. It hence 
takes values in the hermitian 4x4 matrices, and for each (p, x) possesses the two doubly 
degenerate eigenvalues 



h±(p, x) = e<f>(x) ± V \cp - eA(x)) 2 + m 2 c 4 (2.4) 
with associated projection matrices 7r^(p, a;) to the h±(p, cc)-eigenspaces in C 4 . The func- 



tions ( p.4| ) can be identified as the classical Hamiltonians of relativistic systems with pos- 
itive or negative energy, corresponding to particles or antiparticles, respectively. They 
generate the Hamiltonian flows on the phase space R 3 x R 3 , representing the transla- 
tional part of the classical analogue to the quantum dynamics generated by ( [2.1|) . 

It will turn out that the spin part of the relevant classical dynamics is determined by 
the equations 

s±(t) = C ± (& ± (p,x))xs ± (t) (2.5) 

that describe the Thomas precession of a classical normalised spin s 6 S 2 along the tra- 
jectories $±(p, x). Here C±(p,x) contains the electromagnetic fields and potentials as 
known from the relativistic Thomas precession [T5|]. For each classical Hamiltonian h± 
the dynamics of the translational and the spin degrees of freedom can be combined into a 
single (non-Hamiltonian, skew product) flow 

Y&p, x, s) = ($* ± (p, x), 8 ±{t)) , with Yl(p, x, s) = (p, x, s) , (2.6) 
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on the combined phase space R 3 x R 3 x S 2 . On Q^, x S 2 , where = h^}(E) C R 3 x 
an energy shell of h±, the flow Y± leaves the product measure dfj, E that consists of Liouville 
measure on Q E and the normalised area measure on S 2 invariant. If the energy shells are 
compact we will always assume that their Liouville measures are normalised. Then we 
can consider situations in which the flows Y± are ergodic with respect to the probability 
measures dfx E . 

There exists an extensive calculus for Weyl operators of the form ( [2.3|) when their 
symbols obey certain restrictions, see e.g. |L3|, |7|]. If, in the present case, |A(ac)| is bounded 
by some power \x\ , one can introduce the order function 

M(p, x) := v / (cp-eA(a;)) 2 + m 2 c 4 . (2.7) 
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It defines the symbol class S(M), which consists of all smooth, matrix valued functions 
B(p, x) that satisfy the estimates 

\\d$d%B(p,x)\\M < C aP M(p,x) , (2.8) 

where || • 1 1 4x4 is some matrix norm. The symbol Hr>(p,x) of the Dirac-Hamiltonian then 
is in S(M), if the potentials 4>(x) and Ak(x) (k = 1,2,3) are smooth and <f>{x) as well as 
all derivatives of 4>{x) and Ak(x) are bounded. In this case the estimate 

|| (H D (p, x) + i)- 1 || 4x4 < C M(p, x)- 1 (2.9) 

holds and implies that for small enough h the Dirac-Hamiltonian Hd is essentially self- 
adjoint on the domain Cg°(R 3 ) ® C 4 . It therefore generates a unitary time evolution. 

3 Semiclassical projections 

On the level of Weyl symbols there exist the two projection matrices ir^(p, x) onto the two- 
dimensional eigenspaces of H D {p, x) in C 4 . These lead to the two classical flows introduced 
in the preceeding section that correspond to positive and negative energies, respectively. 
In the quantum system a related division of states into particles and antiparticles would 
also be desirable. Due to the effects of pair creation and annihilation such a separation, 
however, can only possibly be achieved in the semiclassical limit. Weyl quantisation sug- 
gests op^^] as first candidates for the desired quantum projectors, but these (bounded) 
operators in fact are no projectors on the Hilbert space 7i. One therefore proceeds to an 
inductive construction by correcting the symbols order by order in h, starting with 7r^ as 
the lowest order term, 

oo 

ft* := op^fvr^ with n ± (p,x) ~ ^h k ^(p,x) , (3.1) 

k=0 

see ||. The expansion in powers of h has to be understood in the sense of an asymptotic 
series m the symbol class S(l) that consists of all B G C°°(R 3 x R 3 ) <g> C 4x4 with bounded 
derivatives, \\d^d^B(p, £c)||4 X 4 < C a p. The Weyl quantisation of such a symbol then is 
a bounded operator on TC [|fl. The terms in the asymptotic expansion fl3.1| ) are uniquely 
fixed by the requirement that (in operator norm) 

(n ± )* = fi ± , IKn^-rPU = o(h°°) , H^n^U = o(h°°) . (3.2) 

That such almost projection operators exist is guaranteed by: 

Proposition 3.1. Let A^ (k = 1,2,3) and <fi be smooth potentials such that and all 
derivatives of A^ and are bounded. Then there exist two bounded almost projection op- 
erators IT that fulfill (\3.3j) and that are Weyl quantisations of symbols tt G S(l) with 
asymptotic expansions ( \3.J\ ). Moreover, the almost projectors provide a semiclassical res- 
olution of unity on 7i, in the sense that ||IT + + Il~ — 1-^11 — 0(h°°). 
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A similar result can be found in fT2" |. Following |]TT| one can now construct genuine 
projection operators P ± from U ± through the Riesz formula, 




These projectors are themselves Weyl operators with symbols in S(l) and are semiclassically 
close to the almost projectors, 1 1 — 1 1 = 0(h°°). However, the operators P^ in general 
do not commute with the Dirac-Hamiltonian. They only almost commute with Hd in the 
same way as the operators Il ± . Nevertheless, one can relate the semiclassical projectors 
P^ to certain spectral projectors of Hd- Depending on the behaviour of the potentials at 
infinity there exist constants E± such that the spectrum of Hd inside (£L, E + ) is discrete; 
e.g., if the potentials and their derivatives vanish at infinity, E± = ±mc 2 . For E G (-E-, E + ) 
we then introduce the spectral projectors P^ ec and P^ pec to the spectral intervals (—00, E) 
and (E,+oo), respectively. 

Proposition 3.2. If E £ (E_,E + ) is neither in the spectrum of Hp, nor an accumula- 
tion point thereof, and if the same conditions as in Proposition [O] hold, the semiclassical 
projectors are close to the spectral projectors, || _P ± — -P S p ec || = 0(h°°). 

The projectors that almost commute with the Dirac-Hamiltonian allow to introduce the 
subspaces := P ± H that in the above semiclassical sense can be viewed as particle and 
antiparticle spaces. The time evolution generated by Hd leaves these subspaces invariant 
up to the semiclassically long time scale T(KJ = 0{h~°°). 

A further consequence of the projectors P ± to almost commute with the Dirac-Hamil- 
tonian is that the projections of eigenvectors ip n G TC of H D , Hd^u = E n ip n , are quasimodes 
with discrepancies of 0(h°°). This means that P ± ^ n G H± are almost eigenvectors of Hd 
with errors whose 7i-norms are oiO(h°°): IKHd-E^P^J =0(H CO ). Appropriate phase 
space lifts of the normalised quasimodes 0j := P ± ip n /\\P ± ipn\\ can now be studied in the 
semiclassical limit and can be related to the dynamical behaviour of the respective classical 
flows K£. 

4 Semiclassical time evolution 

We now consider the time evolution of observables generated by Hd- For convenience we 
restrict attention to bounded Weyl operators B = oj> [B], with symbols from the class 
S(l) that possess an asymptotic expansion in integer powers of h, compare (|3.1|) . We call 
such operators semiclassical observables. The semiclassical projectors P^ are of this type 
and provide a natural separation of an observable B into diagonal and off-diagonal blocks, 
B = B d + B od + 0(h°°), with 

B d := P + BP + + P~BP- and B od := P + BP~ + P~ BP + . (4.1) 
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One expects that the diagonal blocks will be propagated semiclassically by the classical 
dynamics associated with the eigenvalue functions h + and 7t_, respectively. For the off- 
diagonal blocks it is not so obvious how a semiclassical propagation works. In fact, under 
the quantum time evolution the off-diagonal blocks will cease to be semiclassical observ- 
ables. A related discussion can be found in |J. 

For a precise statement we require in addition to the assumptions on the potentials 
made previously that their first and all higher derivatives are bounded: 

Proposition 4.1. Let B = op w [B] be a semiclassical observable. Then for t > its time 
evolution B(t) generated by the Dirac-Hamiltonian is a semiclassical observable, B(t) = 
op w [B(t)} with symbol B {t) G S(l) and asymptotic expansion of the type (pO[], if and only 
ifB od = 0{K°°). 

Hence the propagation of the diagonal blocks can be analysed semiclassically. To leading 
order this will happen in terms of the Hamiltonian flows $± for the translational degrees 
of freedom and the spin-transport matrices D± 6 SU(2) that follow from the equation 

D±(p,x,t) + ±C±(&±(p,x)) -crD ± (p,x,t) = , D ± \ t=0 = 1 2 . (4.2) 

With suitable isometries V±(p,x) : C 2 — > 7T^(p, a;)C 4 these spin transport matrices can 
be converted to d±(p,x,t) e U(4) via D ± (p,x,t) = V±(<l4(p, x))d±(p, x, t)V±(p, x). By 
further introducing the notation B^ := P ± BP ± = op w [B ± ] we can now state the relevant 
Egorov theorem: 

Theorem 4.2. Under the conditions specified above the quantum time evolution of the 
diagonal part Bd of a semiclassical observable B = op w [B] is again a semiclassical observ- 
able, B d {t) = op w [B d (t)], with symbol B d {t) e S(l) and asymptotic expansion 

oo 

B d (t)(p,x) ~ ^2h k B d (t) k (p,x) . (4.3) 

The h-independent, leading term is completely determined by the Hamiltonian flows 
generated by the eigenvalue functions h±, and by the unitary spin-transport matrices d±, 

B d {t) (p,x)= d-\ P ,x,t)B^(^ u (p,x))d u (p,x,t) . (4.4) 

^e{+,-} 

The two types of dynamics that enter on the right-hand side of ( f4.4|) can be combined 
into skew product flows on R 3 x R 3 x SU(2) given by 

Y±(p,x,9) ■= {&±(p,x),D±(p,x,t)g) . (4.5) 

The double covering map R : SU(2) — > SO(3) now allows to relate these flows to the 
genuinely classical skew product flows (|2.6|) , if one sets s±(t) = R(D±(p, x, t))s. Moreover, 
the dynamical properties of the flows Y± and Y± are closely related. E.g., Y± is ergodic 
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with respect to the product measure d/j, E on £l E x S 2 , if and only if Y± is ergodic on 
C!, E x SU(2) with respect to the product of Liouville and (normalised) Haar measure. Thus, 
to leading semiclassical order the time evolution of block-diagonal observables, generated by 
the Dirac-Hamiltonian Hp, can be completely described by the two classical flows Y± that 
combine the translational and the spin degrees of freedom of particles and antiparticles, 
respectively. 



5 Quantum ergodicity 

In quantum systems without spin quantum ergodicity means that phase space lifts of 
almost all eigenf unctions of the quantum Hamiltonian semiclassically converge to Liouville 
measure on an appropriate energy shell, if the classical Hamiltonian flow on this energy 
shell is ergodic. Apart from an Egorov theorem, a proof of this statement requires a 
(Szego-type) limit theorem for averaged phase space lifts of eigenfunctions. 

In the case of a Dirac-Hamiltonian one first has to ensure the very existence of eigen- 
spinors ip n G Ti. To this end we require that there exists an energy E such that all energy 
shells fi^, are compact when E' varies in [E — e,E + e] with some e > 0. These E' shall 
moreover be no critical values of the eigenvalue functions (classical Hamiltonians) h±(p, x). 
In addition, at least one of the energy shells Q E shall be non-empty. Then, for sufficiently 
small h, the spectrum of H D is discrete in the interval I(E, K) := [E — Hlj, E + tiw], com- 
prising of Nj > eigenvalues. 

The desired limit theorem is concerned with averages of the expectation values of a 
semiclassical observable B = op w [B] in normalised eigenstates ip n of H D with eigenvalues 
E n G I(E, h). For this to hold we require that the periodic orbits of the Hamiltonian flows 
<l>± on il E are of Liouville measure zero; e.g., this condition is fulfilled if the flows are 
ergodic. 

Proposition 5.1. Under the conditions imposed above on the Dirac-Hamiltonian the num- 
ber Nj of eigenvalues in the interval I(E, h) grows semiclassically according to 

N - = T° X Tny a ~ E (l + 0{h)) - (5 ' 1) 

Moreover, the Szego-type limit formula, 

lim J_ v <* 6t) i ^+.-^'^W 1 ^ (52) 

K/V, 2^ 2 vol!i++vol!ii ' (5 ' 2 ' 

holds for any semiclassical observable B = op w [B}. 

E v 

Here (tTqB 7Tq) ' denotes an average over VL V E with respect to Liouville measure. On 
average, therefore, expectation values of observables in eigenstates of the Hamiltonian 
semiclassically converge to a weighted average of the 'classical' projections -KqBq-Kq of the 



7 



'classical' observable Bq(p,x); the weights being determined by the relative volumes of 
the respective energy shells. The factor ~ accounts for the two dimensions of the spin 
space (i.e. the range of 7Tq in C 4 ). An important consequence of the limit formula ( |5.2|) is 
that only the diagonal part Bd of an observable gives a non-vanishing contribution to the 
semiclassical average. 

To prove quantum ergodicity now requires to combine Szego-type limits with the 
Egorov-Theorem [O . The latter, however, is only concerned with block- diagonal observ- 



ables. We therefore consider only observables of the type P V BP V . For these Proposition |5Tl 
relates expectation values of B in the projected eigenspinors P u ip n to the classical weighted 
averages of 71qBqttq. At least some of these projected eigenspinors can possibly vanish as 



h — > 0. But Proposition |5.1| allows to conclude that a positive fraction of them retains a 
positive norm in the semiclassical limit; hence these can safely be normalised. As discussed 
in section ^| the projected (and normalised) eigenspinors in general only yield quasimodes 
for the Dirac-Hamiltonian with discrepancies 0(h°°). Thus, quantum ergodicity in the 
present context is concerned with quasimodes rather than with actual eigenstates. The 
reason for this lies in the fact that only after the projection can one associate a definite 
classical dynamics to eigenspinors. 

Theorem 5.2. Suppose that all the conditions hold that have previously been imposed on 
the Dirac-Hamiltonian Hp, as well as that the skew product flow Y* (v G {+, — } fixed) 
defined in eq. ( $.(\ ) is ergodic on Vt u E x S 2 . Then in every sequence {</>^} n eN of normalised 
projected eigenspinors of Hp, with HP^^H > S (8 > small enough and fixed), there exists 
a subsequence {0,^.}jeN of density one, i.e., 

ton * {r ' g^A^l = 1 , (5.3) 
«-» #{n; ||P"<J„|| > S, E„eI(E,K)} 

such that for every semiclassical observable B = op^f-B], 

lim^.,^.) = ±tr(7r^o<)^ • (5.4) 

The subsequence {4>n 3 }j& can ^ e chosen to be independent of the observable B. 

The statement of this theorem can be made more transparent, if one chooses an explicit 
phase space representation of the quasimodes 0*. E.g., if one introduces the matrix valued 
Wigner transform 

W[if>] (p, x) := J e-*P y ip( x - \y) <g> ip{x + \y) dy (5.5) 
of ijj G Tl, the expectation value of B in a state $ v n reads 

(C ^ )= (2^F j j H W Wn]{P^){<B^){p,x)) dpdx{l + 0(h)) . (5.6) 



8 



The conclusion ( |5.4j ) of Theorem |5.2| can therefore be rephrased in terms of the Wigner 
transform (|5.5| ) as 

E5 (dw w ^"^ = hda?^') ~ E) 12 ■ (57) 

Here the limit along the density-one subsequence {4>n-}jeN nas to be understood in a 
weak sense, namely after integration with HqBottq. Thus, the Wigner transforms of the 
quasimodes converge to a uniform distribution on the energy shell Vt v E . In the matrix 
aspect, which represents spin, the result is an equidistribution of 'spin up' and 'spin down'. 
With an appropriate Weyl calculus for spin the latter can also be rephrased in terms of an 
equidistribution on the unit sphere S 2 , see 0, |2j. 



In Theorem [5]2] it was supposed that one of the two skew product flows Y* is ergodic; 
no assumption was made about the complementary classical system. If, however, the other 
energy shell at E is empty, the norms HP^^H must converge to one. Then the statement 
of the theorem applies to a density-one subsequence {ip n }jeN of all eigenspinors ip n with 
E n G I(E, h) in the form 

-Em 



\M nj ,P v BP^ nj ) = \tT{K»B Q K») ' . (5.8) 

In this case quantum ergodicity is hence concerned with the eigenspinors themselves. For 
a Dirac-Hamiltonian this situation is not untypical, since both the particle and the an- 
tiparticle energy shell to be non-empty at the same energy requires very strong potentials, 
with magnitudes comparable to the rest energy mc 2 . But then also the description of a 
relativistic quantum system in a single-particle framework begins to become questionable. 
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